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Abstract 

The aim of this paper is to provide a method for exphcit computation of the Bures 
metric over the space of A''-level quantum system, based on the coset parametrization of 
density matrices. 
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1 Introduction 

Recent developments in the emerging field of quantum information theory have received a 
great deal of attention in investigation of the properties of the set of density matrices of 
an A^-level quantum system. In view of such considerable interest a lot of work has been 
devoted to describe and parameterize density matrices. Any density matrix of an A^-level 
quantum system can be expanded in terms of N'^ — 1 orthogonal generators of SU(N) [1], 
which is a generalization of the Bloch or coherence vector representation for two-level systems. 
Boya et al [2] have shown that the mixed state density matrices for iV-level quantum systems 
can be parameterized in terms of squared components of an [N — l)-sphere and unitary 
matrices. By using the Euler angle parametrization of SU{?>) group [3], Byrd and Slater [4] 
have presented a parametrization for density matrices of three-level systems. An Euler angle- 
based parametrization for density matrices of four- level systems is also introduced in [5] , and 
has been used by Tilma and Sudarshan [6] in order to study the entanglement properties of the 
two-qubit system. A generalized Euler angle parametrization for SU{N) and U{N) groups 
is given by Tilma and Sudarshan [7, 8]. In a comprehensive analysis [9], Zyczkowski and 
Slomczynski analyzed the geometrical properties of the set of mixed states for an arbitrary 
iV-level system and classified the space of density matrices. Di^a [10] has provided an explicit 
parametrization for general iV-dimensional Hermitian operators that may be considered either 
as Hamiltonian or density matrices. The parametrization is based on the factorization of 
N X N unitary matrices [11]. A parametrization useful to study the entanglement properties 
of two-qubit density matrices is also introduced in [12], in which authors have shown that 
the space of two-qubit density matrices can be characterized with 12-dimensional (as real 
manifold) space of complex orthogonal group 5'0(4, C) together with four positive Wootters's 
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numbers [13], where, of course, the normahzation condition reduces the number of parameters 
to 15. 

Efforts have been also made to study the geometry of density matrices. In recent years, 
the Riemannian Bures metric [14] has become an interesting subject for the understanding 
of the geometry of quantum state space. It is the quantum analog of Fisher information in 
classical statistics, i.e. in the subspace of diagonal matrices it induces the statistical distance 
[15]. The Bures measure is monotone in the sense that it does not increase under the action 
of completely positive, trace preserving maps [16]. It is, indeed, minimal among all monotone 
metrics and its extension to pure state is exactly the Fubini-Study metric [16]. The Bures 
distance between any two mixc^d states pi and p2 is a function of their fidelity F{pi, P2) [17, 18] 

2 

(1) 

Fidelity allows one to characterize the closeness of the pair of mixed states pi and p2, so, it 
is an important concept in quantum mechanics, quantum optics and quantum information 
theory. An explicit formula for the infinitesimal Bures distance between p and p + dp was 
found by Hiibner [19] 

d.(p,P + dpf = -5:4^^, (2) 

where \j and jA^), (j = 1, 2, • • • , N) represent eigenvalues and eigenvectors of p, respectively. 
Dittmann has derived several explicit formulas, that do not require any diagonalization pro- 
cedure, for Bures metric on the manifold of finite-dimensional nonsingular density matrices 
[20, 21]. 

The probability measure induced by the Bures metric in the space of mixed quantum 
states has been defined by Hall [22] . The question of how many entangled or separable states 
are there in the set of all quantum states is considered by Zyczkowski et al in Refs. [23, 24]. 
Sommers and Zyczkowski [25] have computed the Bures volume of the (N'^ — l)-dimensional 
convex set and the {N^ — 2)-dimcnsional hyperarea of the density matrices of an A^-levcl 
quantum system. In a considerable work, Slater investigated the use of the volume elements 
of the Bures metric as a natural measure over the {N'^ — l)-dimensional convex set of A^-level 
density matrices to determine or estimate the volume of separable states of the the pairs of 
qubit-qubit [26, 27] and qubit-qutrit [28, 29]. 

By using the Dittmann formula [20] and the Euler-angle parametrization [4], Slater has 
computed the Bures metric for the eight-dimensional state space of the three-level quantum 
systems [30] . In a similar work, but instead in terms of the coset space parametrization, we 
have very recently given an explicit expression for the Bures metric of the space of three-level 
quantum system [31]. The coset parametrization provides a geometrical description of the 
set of density matrices and, as well as the Euler angle parametrization does, eliminates any 
overparametrization of the density matrix [10] . 

The aim of this paper is to provide a method for computation of the Bures metric in 
AT- level quantum systems, based on the Hiibner formula and the coset parametrization of 
density matric;cs. We also use the possibility of factorizing each coset component in terms 
of a diagonal phase matrix and an orthogonal matrix [11, 32]. The paper, therefore, can 
be regarded as a further development in the explicit computation of the Bures metric. We 
show that in the canonical coset parametrization, the Bures metric matrix is divided into two 
blocks, an {N — l)-dimensional diagonal matrix corresponding to the eigenvalues coordinates, 



dB{pi,P2) = j2-2^F{pi,p2), F{pi,p2) = 
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and an N{N — l)-dimensional matrix corresponding to the coset coordinates. It therefore 
provides a factorization of the Bures measure on the space of density matrices as the product 
of the measure on the space of eigenvalues and the truncated Haar measure on the space of 
unitary matrices. It is shown that the coset parametrization gives a compact expression for 
all metric elements. The analytical expression for Bures metric can be used in computing the 
Bures volumes of quantum states as well as to study the problem of what proportion of the 
convex set of the bipartite systems is separable. The results also enable the calculation of the 
minimum Bures distance of a given density matrix from the convex set of separable states, in 
order to quantify entanglement of the state [33, 34]. 

The paper is organized as follows: In section 2, the coset space parametrization of an 
A''-level density matrix is introduced. Based on the parametrization, we provide in section 3 
a formula for computation of the Bures metric, explicitly. The paper is concluded in section 
4 with a brief conclusion. 

2 Canonical coset parametrization of density matrices 

The state space of an A^— level quantum system is identified with the set of all TV x Hermitian 
positive semidefinite complex matrices of trace unity, and comprises an {N^ — l)-dimensional 
convex set A4n- The total number of independent variables needed to parameterize a density 
matrix p is equal to iV^ — 1, provided no degeneracy occurs. Let us denote the set of all diagonal 
density matrices of an AT- level quantum system by V^- An arbitrary element p^^^ € "Dat can 
be written as 

N 

p(^) = diag{Ai,A2,--- ,Ajv}, 0<Ai<l, ^Ai = l, (3) 

1=1 

which simply denotes an {N — l)-dimensional simplex Sn-i for the set of all diagonal density 
matrices. 

A generic density matrix p G jv in an arbitrary basis can be obtained as the orbit of 
points p^^^ S X>jv under the action of the unitary group U{N) as 

p=Up^'^^Ul (4) 

Let H he a maximum stability subgroup, i.e. a subgroup of U (N) that consists of all the 
group elements h that will leave the diagonal state p^^^ invariant, 

M''^/it = p(^), heH, p(^)eCjv, (5) 

that is, H contains all elements of U{N) that commute with p^^\ For every element U G 
U{N), there is a unique decomposition of U into a product of two group elements, one in H 
and the other in the quotient G/H [35], i.e. 

U = nh, UeU{N), h€H, n&U{N)/H. (6) 

Therefore the action of an arbitrary group element U e U{N) on the point D G is given 

by 

p = up^'^^u'' = mp^'^^h^si^ = np(^)nt. (7) 

Since Djv consists of points with different degree of degeneracy, the maximum stability sub- 
group will differ for different p^^^ G Vn [9]. Let rrii denotes degree of degeneracy of eigenvalue 
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Ai of matrix p^^^ . It follows from this kind of spectrum that p^^^ remains invariant under the 
action of arbitrary unitary transformation performed in each of the mj-dimensional eigensub- 

spaces. Therefore H = (7 (mi) (g) U{m2) (X" • • • U{mk) is maximum stability subgroup for p'^^\ 
and the quotient space U{n)/H is a complex flag manifold 

^= U{m,)^Uim,)^...Uim,y + -^2 + ■ ■ ■ + m, = N. (8) 

Two special kinds of degeneracy of the spectrum of p^^^ are as follows: i) Let p^^^ represents 
the maximally mixed state = diag{-^, jfr'':jf}- In this case the stability subgroup H 
is U{N), and the orbit of point p* is only one point, i.e. p = p^,. ii) On the other hand if 
the spectrum of p^^^ is non-degenerate, then the stability subgroup is n-dimensional torus 
= J7(l)®^, and the orbit of the point p(^) is 

p = J2p(^)J2t, J2 g U{N)/T^. (9) 

Since the maximal torus is itself a subgroup of all maximum stability subgroups, therefore 

the orbit of points p*^^) G under the action of quotient U{N)/T^ generates all points 
of the space Mn- The diagonal matrix p^^^ is defined up to a permutation of its entries 
and, one can divide the simplex <Sjv-i into n! identical simplexes and take any of them. Each 
part identifies points of Sn-i which have the same coordinates, but with different ordering, 
and can be considered as the homomorphic image of simplex Sn-i relative to the discrete 
permutation group Pn, i.e. Sn-i/Pn- Therefore the points of A4n can be characterized as 
the orbit of diagonal matrices D <E Sn-i/Pn under the action of quotient ft G U{N)/T^ . 
Further insight into the space of density matrices can be obtained by writing the quotient 
e U{N)/T^ as the product of A'' - 1 components as ([35], page 401) 

n = Q{N;N)^{N-l;N) . . . j^(2;JV) _ (^q) 

where 

A typical coset representative f2(™;^) can be written as 

= ( SU{m)/U{m-l) \ O \ ^ ^^^^ 

where O, and In-tu represent, respectively, the mx (N—m) zero matrix, its transpose and 
the (N—m) x {N—m) identity matrix. The 2(m— l)-dimensional coset space SU{m) /U{m— 1) 
has the following m x m matrix representation ([35], page 351) 



SU{m)/U{m-l) = 



cos -v/S(™) [BM]t 

anV[B("')]tB('^) rp(mm 
^[B(™)]tB(-) ^ J 



iy[B('")]tB('") 



cos v1b(™)]Tbm 



(13) 



where B^™) represents an (to — 1) x 1 complex matrix and [i?(™)]^ is its adjoint. 

Now by parameterizing the column matrix as (7i™''e*^i ' , 72'"''e*^2 ' , • • • , 7^\e'^"-i 
for m = 2,3, • • • , A'', where j^"^^ and 4^"*^ ai"© real numbers-^, the component Ct^^'^^ can be 

^The correspondence between general notation of this paper and that in Ref. [31], where the N = 3 case 
is explicitly computed, is as 7^^^ = a, = 4>, = /3i,2 and = V'1,2- 
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factorized as 

^{m;N) ^^im;N)j^(m;N)j^(m;N)f m = 2,3,---,N, (14) 

where X^™'^) is a diagonal N x N phase matrix with xj^'^'' = Skicxjp{i^^^} and ^j"*^ = 
for i > m, and R("^'^^ is an N x N orthogonal matrix with the following nonzero elements 

= Sij + 7^'7i™^ (cos i^"^ - 1) for l<i,j<m-l 
lt:''^ = -B^''^=t^^l^ for l<^<m-l 

i?(;r^^) = cos7('") 

^(m;JV) ^ ^ m + 1 < i < AT (15) 

where we have defined 7!*"^ = 7|'"V7^™^ ^-^^d 7^'"' = \j {ll"^^ Y ■ As we will see later 
the important ingredient of our approach in computing the Bures metric is the possibility of 
writing the factorization (14). 



3 Bures metric 

In this section wc shall attempt to develop a method of computing the Bures metric of an 
arbitrary density matrix of an A^-level quantum system. We will use the canonical coset 
parametrization of the density matrices introduced in the last section. 

Let (0 be a generic density matrix of an TV-level quantum system, with eigenvalues Aj 
and corresponding eigenvectors |Aj), {j = 1,2, ••• ,N). In the light of p = flp^^^fl\ with 
p(^) = diag{Ai, A2, • • • , Ajv} as the diagonal matrix of p eigenvalues, the p eigenvectors can be 
written in terms of p*^^^ eigenvectors as jA^) = fl\i). Therefore invoking the Hiibner formula 
(2), we can write the infinitesimal Bures distance between p and p + dp as 

^,,2 1 ^ |((i|»t) d(^lp(^)^|t)(n|j))|^ 

dB(p,p+dp) = - y: xrrx ' ^'^^ 

i,j=l ■> 



which takes the form (remember that CI is unitary and therefore dSl^ = —Cl^ dfl fi^) 

N 

2 A, + Xj 

By using the equations p^'°^|i) = and = we get 

^ 'dA^)^ ^ 
4Ai 



da(.■.+dp)-=iI: l"|^'''°'|^•>^'^f:^"■^'°'||^'|^ (17) 



ds(p, p + dp)2 = ^ i-^ + ^ A,,|(ntdn),,f , (18) 



i<3 



where we have defined as 
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g(D) 











Equation (18) shows that the infinitesimal Bures distance is divided into two infinitesimal 
distances corresponding to the eigenvalues coordinates and coset coordinates. Therefore the 
Bures metric matrix becomes block diagonal as 

(20) 

where g^^^ is a part of the Bures metric corresponding to the diagonal density matrix p^^\ 
and g^'-'^ is the contribution of the coset coordinates in the Bures metric. 

In what follows our goal is to calculate the matrix elements of g^^^ and g^'-"^ . In order to 
calculate g^^\ we first note that the N eigenvalues can be parameterized explicitly in terms 
of A'' — 1 independent parameters 9k {k = 1, 2, • • • ,N — 1) as 

A/j = sin^ ^1 sin^ ^2 • • • sin^ Ok-i cos^ 6k for k = 1,2, - ■ ■ ,N — 1 

(21) 

\n = sin^ 9\ sin^ 62 - ■■ sin^ On-i- 

By using the above coordinates for the eigenvalues of p, we can write the infinitesimal Bures 
distance over the eigenvalues coordinates as 

N /, , n2 N-l 

where it can be easily seen that the metric gj.^^ is diagonal with elements 



-■^^^ -sin26'isin2 6'2---sin2 6'fc-i, for A; = 2, 3, • • • , iV - 1. (23) 



It is worth noting that the metric g*^^^ is independent of 9n-i- 

Now in order to calculate g^^\ we define {<^m = 1, 2, • • • 1 2(m — 1)) as the 2(m — 1) 
real parameters of the coset component fi^™'^) (m = 2, • • • , A'') such that 

for a„ = 1, • • • , m - 1, 

(24) 

for a„ = m,- • • ,2(m- 1), 
Then the infinitesimal Bures distance over the coset coordinates can be written as 

= EL2El":,-:li5i3r^dx«dx« (25) 

, ,p2(m-l)v^2(m'-l) (m,m';JV), (m), (m') 

Now invoking the decomposition 

f^^^iN;N)^iN-l;N),,,^i2;N)^ (26) 
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we can write 



N 2 N f ^ 

\{nUn)ijf = + 2 ^ Re | (^K^rn'-.N)^ 1 (27) 

m=2 m<m' 



where we have defined 



ith 



W{m;N) ^ Q(m-l;JV) . . . q(3;N)q{2;N) 



m;N) 



(28) 



(29) 



and = 1. The first term of the r.h.s. of equation (27) shows sum over all pure 

contribution of each coset component in the Bures distance, and can be identified with 



2(m-l) 



(30) 



The second term, however, shows the sum over mixed contribution of all pairs of coset com- 
ponents in the Bures distance and can be used to write 

n ( s 2(m-l) 2(m'-l) 

EA,,Re|(i^(-^))^.(i^(-'^^))* 1= E i'if^dxt^dx^?. (31) 



Therefore ^f^'-^-' is defined by the following symmetric matrix 



/ 9 



(2,2;JV) 



V" 



9 



(2,3;Af) 



J3,3;Af) 



9 



(2,N-N) 



9 



(32) 



Now the object is to calculate the matrix elements of g^*^^ . which can be achieved if we can find 
an expression for (if ^™'^-')^^. of equation (30), (31) in terms of the coset coordinates X^m {or 

equivalently in terms of 7!™'' and ^i™'). Therefore we have to expand r2(™'^)^df2'-™'^^ of 
equation (28) in terms of d7(™) and d^(™). This can be achieved by using the factorization 

Q(m;JV) ^ ^(m;iV)^(m;iV)^(m;iV)t ^j^^^.^ 

Q(m;JV)t^Q(TO;iV) _ j^(m;JV)^(m;JV)'^j^(m;iV)t^^j^(m;iV)^^(m;iV)j^(m;JV)t 

- (dX(™;^))X(™;^)^ (33) 

+ Jf(m;Ar)^(m;JV)'^^jj^(TO;iV)-jj^(TO;iV)t_ 

By knowing that X^"^'^^ is a unitary diagonal matrix and is an orthogonal matrix, 

the first two terms can be calculated easily by using the relation 



dxg"^^)=ie'€^'5.,dC^\ 



(34) 
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On the other hand to calculate the third term, the equation (15) may be used to show that 



m — 1 



di?lr"'=ErS"^d7r, (35) 



where for a fixed fc (= 1, 2, • • • , m — 1), the ^\j!if^^ denote the matrix elements of an N x N 
matrix with nonzero elements 

^ (m) ^ (m) ^ (m) . (m) r i ^ ■ ■ ^ i 

~li % % sm7*- for l<«,j<m— 1 



p(m;JV) _ _p(m;JV) 



(e ,(m),(m)\ . (m) / (m) 

^ik-li Ik ) sin7^ '/t ' 

_^ . (m) 1 < i < TO - 1 

rS = -^f^«in7(-) (36) 
Using equations (34) and (35) in equation (33), we get 

f0('":^)^df2(™'^)) = expO(el"''^^ - d"'"^^)} 



f(m) 



r=l 

where we have defined 



(37) 



where ^-E^J^'^^j and (^E^J^'^^^ have been defined by 

m 

=E^lr^"^r{;^-) (38) 

(e^p^a =4r'^)4r^)+5,,^,. (39) 

Therefore for if ('"s^) of equation (28), we find the following matrix elements 

m— 1 J, 



= E K^^'T • (41) 

and 

= E (^r^"T (^^"^"0. ' (42) 
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Finally, putting everything together, we find the following formula for the matrix elements of 
the metric g^*^^ 



N 



^,.,..,-^^..Re{(Alr"'),w(<'"');-.}- («) 

i<j 

for r = 1, 2, • • • , m — 1, s = 1, 2, • • • , m' — 1 and m, m' = 1, 2, ■ • • , A''. Since the eigenvalues 

coordinates arc just included in the A, j terms of equations (43)- (46), therefore it is clear that 
all matrix elements of the metric g^'~'^ are simply sum of the products of two independent 
functions, the eigenvalues coordinates and the coset coordinates. It should be noted that 
for real density matrices, i.e. ^j-™-* — for i = 1, 2, • • • , m — 1 and m = 1, 2, • • • , A^, we 
have = R^"^'-^\ In this case all quantities are real and equation (43) gives all matrix 

elements of the N[N — l)/2-dimensional metric g^'^h 

4 Conclusion 

We provide a method for explicit computation of the Bures metric in A^-level quantum systems, 
based on the coset parametrization of density matrices. We show that in the canonical coset 
parametrization, the Bures metric matrix is divided into two blocks, an {N — 1) -dimensional 
diagonal matrix corresponding to the eigenvalues coordinates, and an A^(A^ — 1) -dimensional 
matrix corresponding to the coset coordinates. It also provides a factorization of the Bures 
measure on the space of density matrices as the product of the measure on the space of 
eigenvalues and the truncated Haar measure on the space of unitary matrices. 
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